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INTRODUCTION 



Inspired by String/M-Theory brane-physics have blossomed in recent years. Branes are p dimensional 
extended objects embedded in a higher D > p+1 dimensional universe which arise naturally in String 
theory as hypcrplancs where open strings can end, and in Supergravity/M-Threory as solitons [1]. 
Mostly, branes have found applications in two virtually opposite set ups: In the "brane-world" 
scenario [2] the brane dimensions coincide with the three observed dimensions of our universe, while 
the spatial dimensions perpendicular to the brane, which can be both compact or non-compact 
correspond to the "unseen internal dimensions". Contrastingly in the "brane-gas" scenario [3,4] 
the branes wrap around compact internal dimensions, while the directions perpendicular to the 
branes become the observed dimensions of our universe. While the brane-world picture have several 
virtues including being able to address the hierarchy and the cosmological constant problem, the 
brane gas model can explain why the internal manifold remained small as compared to the observed 
dimensions, and a simple counting argument also yields the dimensionality of the observed universe 
to be three! Invoking T-duality brane gas cosmology (BGC) also seems to be able to avoid the big 
bang singularity. Thus, it may be an interesting venture to combine the two scenarios by considering 
branes some (three) of whose dimensions are noncompact (which becomes our observed universe) 
while the others wrap compact internal manifolds, which potentially can be non-trivially curved - 
the isometries of this curved internal manifold can then be associated with gauge fields living on 
the brane. In this paper however, we only restrict ourselves to the usual set up of BGC 1 * but we 
generalize the internal manifold from being a flat tori to a curved space. 

In [4] Brandenberger and Vafa observed that string winding states would tend to prevent expansion 
since the energy of the states increases if the circumference increases. However, strings with opposite 
winding numbers can annihilate if their world volume intersects and this happens most efficiently 
in three dimensions. Thus in a three dimensional subspace the string winding states can annihilate 
letting these dimensions grow, while strings winding other directions will not be able to annihilate 
each other efficiently and ultimately fall out of thermal equilibrium thereby stopping the expansion. 
This stabilization mechanism of the internal manifold by employing 'string gas' was further quantified 
in [5] and has also been generalized to 'brane gas' scenarios [6] (for more details also see [7]) with 
essentially similar results, except that now we have a hierarchy in the sizes of the extra dimensions, 
coming from the contributions of the p-branes, with different p, as in general p-branes annihilate 
most effectively in 2p+l dimensions. Although BGC has been successful in solving several problems 
associated with standard big bang cosmology, there are some unresolved issues as pointed out in 
[8]. For example, it is not clear exactly how the three dimensions that we observe today are chosen 
among all the other dimensions, perhaps through random quantum and thermal fluctuations. In 
that case it is possible that our universe may contain different patches where different directions have 
become large! Also, it is known that branes with different values of p can interact with each other, 
changing their winding numbers, thereby allowing thermal equilibrium to be maintained among the 
brane winding states. A simple way to eliminate these issues is to assume the universe to be a direct 
product of three non-compact dimensions and a compact internal manifold. The BGC in such a 
set up has been extensively studied in [8] when the internal manifold is a tori. Here, we adopt the 
same approach, but generalize the program to nontrivial curved internal manifolds* to study how, 
if at all, the internal curvature may change the cosmological dynamics. Of course, in this set up the 
question of why the universe is a product space with specifically three non-compact dimensions is 
an open issue, and the nice BGC argument concerning the dimensionality of our observed universe 
becomes redundant. We will come back to this at the end. 

For a flat p-brane the p longitudinal (along the brane) directions are flat, or in other words a flat 
p-brane solution preserves the isometries of the p dimensional flat space. These branes can then wrap 



fin the strict sense our study is slightly different from BGC as the dilaton is fixed to be a constant in our 
case as we will soon find out. However, the analysis is related and the value of the dilaton is supposed to 
get stabilized at later times any way. 

*It should be mentioned that BGC with branes wrapping nontrivial cycles (in Calabi-Yau and K3 spaces) 
was first studied in [9] . 
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around flat p-dimcnsional tori. If one now considers brancs wrapping a curved internal manifold, 
say M p then the brane solution should accordingly, preserve the isometries of M p . In section I, we 
derive such p-brane solutions^ which look like de Sitter black holes from the four dimensional point 
of view. This is to be expected as after reduction ordinary p-branes resemble black hole solutions, 
and since our internal background manifold is curved, we are living in an asymptotically de Sitter 
space-time. We here only obtain the uncharged solution but it may be interesting to study other 
charged solutions which will have underlying connections to non-abelian black holes (see [11] for 
a review). In section II, we obtain the back reaction of these branes on the space-time geometry 
by analogy with the flat brane solution [8] and check that the stress energy tensor is consistent. 
As is usual the branes represent delta function sources for the Einstein's equations which can be 
smoothed out by considering a uniform density of brane gas. In section III, we study BGC both 
without and with matter/radiation. We find that the early evolution resembles the picture with flat 
branes [4,6,8]: while the branes wrapping the internal manifold tends to contract them, the internal 
curvature terms only reinforcing this effect, the noncompact dimensions on the other hand are free 
to expand as the branes act as pressureless dust. At late times however the picture changes as the 
curvature term starts dominating over the brane effects. The stabilization of the internal manifold 
that was achieved in the presence of ordinary matter and branes (the expansion of ordinary matter 
being balanced by the contraction of the branes) is lost. Preliminary analysis indicates that the 
internal manifold may also be growing along with the ordinary dimensions. This suggests that we 
may have to invoke other mechanisms to stabilize the internal manifold at late times, like turning 
on the fluxes [14]. We conclude by summarizing and commenting on possible future research. 

I. BRANES WRAPPING CURVED MANIFOLDS 

Supergravity and Brane Ansatz: We start with a generic bosonic sector** of a supergravity 
(SUGRA) action: 

S = -L_ J dx ^ g{ R _ \ dfh ^ - \ £ _L e ^ _ 2A } (1) 

where (f> is the dilaton field and field strengths F ni 's are nj forms, / = 1 . . . M and we have also 
included a cosmological constant term. Although both in the action and in the field strengths one 
can have Chern-Simons-like terms it is known that they can be ignored while considering brane like 
solutions. The field equations that follow from (1) are 

Gfhn ~ lidfhcp^ - ^gfnnid^) 2 ) - ^ e a ' * (njF^ F^ -w - l 9fiA Fl) + A gfh n = (2) 




We now specialize to the uncharged solution, for which we can consistently put the field-strength 
and the dilaton to zero. We are thus left with only (2y'; (3) and (4) being trivially satisfied. 

Now let us look at the ansatz for our p-brane metric. For simplicity we assume the internal 
manifold to be a group manifold G e , but one should be able to generalize the results to other 

D 



§ Previously branes wrapping curved dimensions have been studied using boundary state conformal field 
theory [10], but here we take a supergravity approach. 
**In brane solutions all fermionic fields are typically put to zero. 
t1 We are essentially in Einstein's gravity with a cosmological constant. 
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homogeneous spaces. Since in our picture the p-branc is wrapping the internal manifold, D= p. 

Moreover, the p-brane metric have to possess the isometries of the internal group manifold. As in 



the flat p-brane case, this fixes the D directional components of the metric: 

(9mn{x) \ 

^{x)g% (y) ( 5 ) 
ran / 

We employ symbols '0' and 'no symbol' to denote quantities corresponding to the internal and 
the external manifold respectively, wherever necessary. Also, x and y denote coordinates charting 
the observational universe and the internal manifold respectively. Here g K is the Killing metric** 
possessing the isometries of the group manifold and ^(x) is the radion. If we now further impose 
invariance under time translations and SO(D) rotations among the transverse x directions then we 
obtain 

ds 2 = ~p 2 dt 2 + Pdr 2 + Y 2 r 2 dtf D _ x +$ 2 ds 2 (6) 

where all the functions only depend on the radial coordinate r, the brane being located at r = 0. 
Note that the only way this metric differs from the usual flat brane ansatz is that the metric along 
the brane directions is not flat but corresponds to the Killing metric of the group manifold. At this 
point one could use (6) to solve the field equations (2) but it is easier to perform a dimensional 
reduction which then maps the p-branes to black holes. Since the internal manifold is curved we, in 
fact, expect the p-branes to look like de Sitter black holes [15]. 

Reduction and Brane Solution: To perform the reduction it is convenient to work with the 
viclbcin: 

. / e m a (x) \ „ ( e a m (x) \ 



and the "flat-metric" 



fh& = (8) 



'©9 

a b 



Such a dimensional reduction is known to be consistent and yields 



S = idol - 2X > = TedcJ 1"'^ J^^TjAr + k - V] 



where 



and 



K=-(2D—+D (D -iy-^r) ; V = ~*~ 2 + 2A 




16ttG s J " v 16ttG s 167rGr. + i 



**One can also employ more general "squashed metric" but we will not consider them in this manuscript. 
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To obtain the action in the Einstein frame, we perform the well known conformal re-scalings 

o 

e a a = Ael ft ; A = (9) 



Then we have 



with 



K' = - Di P 1! ' id'*)'' ; V = 2Ae~ 2 A* - ^ e - 2 * (10) 



where we have defined 

*' = 

In future we will drop the 'primes'. 

The action (10) yields field equations for the four dimensional metric and the radion. However, 
we know that for uncharged p-branes the radion is a constant satisfying 

dVW n 



D-2 



The potential V then acts as an effective D + 1 dimensional cosmological constant A given by 

A.^e-^^fjM^ (12) 
8 8 \D-2j 

The effective action now reads 

De Sitter black hole solutions for (13) are known [15]: 

ds 2 = -f3 2 dt 2 + $ 2 dr 2 + rVdQfj.j (14) 

with 

r = l ; /3 = f>- 1 ; $ 2 = ~ — -r = ^ (15) 

1 ip 2A 2 ffr) 

The ADM mass of the black hole T p , arises as an integration constant and is equivalent to the ADM 
mass or tension of the p-branc [12]. For completeness sake, let us write down the full brane solution 

g / _ 3-2 

^ - (^) D1 -fWt 2 + ^ + rW^ + (^) ^ 1 s 2 | (16) 
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II. BRANE BACK REACTION 



Having obtained the brane solution wrapping the extra dimensional manifold, the next task would 
be to compute its back reaction on the background geometry. The cosmological ansatz for the 
background metric looks like 

ds 2 = -p 2 dt 2 + a 2 dx 2 + f 2 dy 2 (17) 

where we take the external three dimensional spatial geometry to be flat. Again it is useful to 
introduce the viclbcins 



( e~ w \ 

e- A 5u™ 

a 

y e- s eg™y 



(18) 



along with the flat metric (8). We use the symbol 'U' for quantities and indices corresponding to the 
external spatial dimensions. Here we have also redefined f3, a and iff as exponentials of W, A and 
S respectively, which are now to be treated as collective coordinates characterizing the geometry of 
both the internal and the external manifold, depending only on time. In [8] the stress energy tensor 
for a flat brane located at x in such a background was already obtained which lends a natural 
generalization to curved branes §§ : 

Too = T p e- DA S(x - x ) 
T uU = 

ab 

T Q& = -T p e- DA 5(x - x )g @e (19) 

ab ab 

Consider now the brane gas scenario with say ni branes located at xi which implies S(x — xq) — > 
J2 X , ni8{x — xi). As is standard practice we now pass on from a discrete to a continuous distribution 
of brane gas 



^ niS(x — xi) — > j dx'n(x')5(x — x') 



where n(x) is the density of brane gas. Assuming a uniform density then gives us 

Too = nT p e- DA 



T u u = 

ab 

T 00 = -nT p e- DA g ee (20) 

a b ab 

This stress-energy tensor differs from the usual brane gas stress-energy tensor (see for example [6]) 
which has an equation of state parameter —p/(D — 1) along all dimensions. As we will see later, 



^Although the result physically makes sense because we expect the branes when wrapped around the 
internal manifold (curved or flat) to behave as delta function sources in the transverse (observable) directions, 
it would be nice to derive it ab initio starting from, for example, the Polyakov brane action. This however 
would require a better understanding of the brane dynamics which may also lead to some corrections to 
(19). 
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(20) corresponds to equation of state parameters -1 along the wrapped compact dimensions and 
for the non-compact unwrapped directions. This apparent discrepancy arises because we do not 
perform an average of brane wrappings over all dimensions as in our case the wrapped dimensions 
are fixed and topologically distinct from the non-compact dimensions. In the conventional BGC one 
has to perform the average as all dimensions are equivalent and then the averaged equation of state 
parameter will indeed be given by 

N 
(-!).£> +0.D D 



D-l D-l 

Before proceeding any further, let us check that the stress energy tensor is indeed consistent, i.e. it 
satisfies the divergence free condition 

V a T a? = (21) 

Now 

Vrriab „ rriab , , 7l rricb , b ryiaC 

a± — ^a 1 + W ca-i + W ca 1 

and the non-zero connection co-efficients for (18) are given by 

w° u=5u"e- w A 

Oc c 

£ 8 „ = sJe~ w S 



'o = e- s J ee (22) 

be be 



Then trivially 



Finally 



V a T ab = 

VtiSO „ rriaO i , ,a rpcO , , ,0 rpac 
a-1 — +W ffi-1 + W ca± 



= e~ w d t {nT p e- DA ) + e- w {DA+ D S)nT p e~ DA - e~ w nT v £ Se~ DA = 



Q 00 00 



+ LO a Qc ,T db +cu b Qa T ad =0 



00 ± "I" 1 " 00 J 

da da 







since w a 00 = for unimodular groups, and compact groups are unimodular. 

da 

III. BRANE GAS COSMOLOGY 

Without Matter: Having obtained the energy momentum tensor for the brane gas, we can now 
proceed to obtain Einstein's equations of motion. In order to do that we first compute the Einstein 
tensor for the geometry (18,8). 



Goo = e 



-2W 



X -D{D - l)i 2 + \d(D -l)S 2 + DDAS+\d\ e 2 ^"^ 
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G u u 

ab 



-d u ue 

ab 



-2W 



(D - 1)1+ D S + l -D{D - l)A 2 + \d{D +1)S 2 



+ D(D- 1)AS - (D - I) AW- DSW+^DX e 2 ^^ 



G © — —g@ee 

ab ab 



-2W 



DA + (B + \d{D + 1)A 2 + \d(D -l)S 2 



+D(D -I) AS - DAW - (B -l)SW + i(B -2) A e 2{w ~ s ^ 



Einstein's field equations then read 



(23) 



(24) 



To analyse these equations in more detail it is convenient to choose the 'canonical gauge'*** [8]: 



W = DA+ D S 



(25) 



After some recombinations one obtains 



(D - 2)1 = (B +l)K 2 nT p e 2 ° s+DA 



(D - 2)S = -(£> - 2)n 2 nT p e 2DS+DA - \ (D — 2 )e 2 ( D ~ 1 1 s+2DA 



(26) 



The third equation is not expected to be linearly independent as usually happens in general relativity 
though it can constraint the initial conditions. These equations look the same as with branes 
wrapping flat tori [8], the only addition being the internal curvature term. However, the curvature 
term has the same sign as that of the brane term and thus it can only enhance the brane gas effect 
of trying to contract the extra dimensions. Unfortunately we couldn't find an analytical solution as 
it was possible in the flat case, but numerical evolution of the equations exemplify the behaviour. 



***This is the gauge where in the action the kinetic terms for the variables S and A are canonical, and 
hence we expect Newton's law type equations which provide a clear picture of the dynamics. One can always 
go back to the more standard W = gauge. 
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FIG. 1. Typical evolution of A(t) and S(t) in Planckian units. 



As expected we find that while the unwrapped dimensions grow, the wrapped dimensions shrink 
(see figure) thus providing an explanation for the differential size of the external and the internal 
dimensions. 

With Matter: Let us now make the picture more realistic by adding matter/radiation: 

rpM _ 

= S uu p 

ab a b 



ab a b 

One can solve for the energy density p using the equations of state 

p = ujp 



(27) 





and the divergence equation 

V s T ao = 

One obtains exactly the same result as for the flat case: 



p = p e 



-D(1+lj)A-D(1+S)S 



(28) 



(29) 
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It is easy to check that other equations in (21) are automatically satisfied. It is interesting to note 

that the brane stress energy tensor is also of the same form as (27,28) with u = 0, S= —1 and 
po = nTp. 

The addition of radiation/matter modifies the equations of motion to 

(D - 2)1 = (B +l)K 2 nT p e 2 ° s+DA + k 2 Po (1 + (£ -l)w- B§) e ^(i-«M+B(i-S)S 

(5 - 2)S = -(£> - 2) K 2 nT p e 2 ° s+DA - A (D - 2 )e 2 ^ s+2DA 

+k 2 Po (1 + (£> - 1) 8 _^) e ^(i-)^+B(i-S)s (3Q) 

It is clear that the evolution would be complex and there could be several regimes depending upon 
the specific values of the equation of state parameters and the values of S and A. In general one 
cannot solve these equations analytically and to comprehensively understand the dynamics one will 
have to resort to numerical simulations. However, let us consider a likely scenario: 

It is desirable to account for an inflationary phase in the beginning. Thus one can assume that 

initially po comes from a cosmological term (or perhaps an effective cosmological term arising from 



slow rolling of a scalar field, with po = Pinfiaton) with uj =lu= — 1. It is easy to see from the evolution 
equations that when both A and S are small enough then the dominant exponents are associated 
with the internal curvature and the brane contributions and hence the evolution would proceed 
similar to the case when no matter was present. A will expand while S will shrink giving us the 
origin of the differential size of A and S. As A increases, there will come a point when cosmological 
term starts to dominate over the brane term. Also note, the internal curvature term always dominate 
over the brane term because S is small while A is large. Thus the phase of differential 'growth' will 
soon give way to an inflationary phase where the size of the internal manifold remains approximately 
constant while that of the external manifold increases exponentially (this is nothing but the familiar 
de Sitter vacuum). After the end of inflation, presumably the universe becomes radiation and later 
matter dominated: 

Pinflaton * /^radiation * Anatter 

At this point one may wonder as to what kind of equation of state parameter does matter and 
radiation obey along the extra dimensions? Although the answer is not clear it seems natural that 
since the internal manifold is now very small the wave- function of the matter/radiation particles will 
wrap around the extra-dimensional manifold and thus effectively behave in much the same way as 

branes do, or in other words have 8= — 1. Since by this time we expect most of the brancs to have 
annihilated the evolution would be governed by a competition between the internal curvature and 
the radiation/matter terms, unlike between brane and radiation/matter as was discussed in [8] for 
the flat case. 

One can now try to find exact solutions (assuming T p = 0) of (30) by matching the exponents 
[13,8]. We find it convenient to directly work in the W = gauge for this purpose, where the 
independent equations look like 

l -D{D - l)A 2 + \f){D -1)5 2 + D D AS + \ DX e~ 2S = n 2 poe -^(i+-)A-B(i+2)S 

and 

DA + (B -1)5 + ±D(D + l)i 2 + \d(D -l)S 2 
+D(f) -1)AS + 1(B -2) A e- 2S = -n 2 S ^-^1+^^(1+2)5 (31) 
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To match the exponents we make the following ansatz 

A(t) = ln{a ) + ailn(t) ; S(t) = ln(s ) + s x ln{t) (32) 
With (32) since A, S, A, S, all go as ~ t~ 2 we impose 

-£>(1 + w)oi- B (1+ S)si = -2 

and 

-2si = -2 



si = 1 and oi = — — — r (33) 

for £> = 3 and §= — 1. For both radiation (w = 1/3) and matter (uj = 0) this gives the usual 
scaling laws for the size of the external universe, e A <~ t 1 / 2 and e A <~ £ 2 / 3 respectively. Matching the 
exponents however still leaves us to satisfy two equations with two unknowns (a , s ). In the flat 
case with branes and matter/radiation, one could indeed find such a solution set, but unfortunately 
it fails in our case (although solutions exist for u — > 0~). However, as remarked in [8] (33) can still 
be taken to indicate a general behaviour of the solutions and we here take a similar view point. 
We also observe that unlike in the flat case where an interplay between the brane contraction and 
expansion due to pressureless dust could stabilize the internal manifold, the internal manifold here 
grows linearly with time! This is because the interplay of curvature terms and matter gives us a 
"potential hill" rather than a valley To exemplify this consider the evolution of S in (30) with 
T p = and in the presence of a pure cosmological constant^. In the equation of motion of for S 
the A dependence factorizes and the evolution is governed by an effective potential of the form 

V ell (S) ~ C cnrv e 2 ^ s - C cos/m&t e 2 ° s (34) 

Clearly there is an unstable potential hill and because of the cosmological/matter term S tends to 
grow once to the right hand side of the hill. This suggests that one perhaps needs to incorporate 
other forms of stabilization mechanism, for example by turning on the flux fields [14], to be able to 
stabilize the internal manifold when it is curved. 

SUMMARY AND FUTURE RESEARCH 

In this paper we have tried to describe a framework to study BGC with branes wrapping curved 
group manifolds rather than flat tori which has been extensively studied. We find that although 
at earlier times the dynamics is similar to the flat case, the branes now reinforced by the internal 
curvature effects again tries to contract the wrapped compact dimensions, the later dynamics is 
different as the internal curvature terms become important. In particular, there seems to be a 
regime when the competition between the internal curvature terms and the usual matter-radiation 
contributions allow the external universe to scale in the usual way <~ t 1 / 2 for radiation and <~ i 2 / 3 for 
matter at least approximately. However the analysis carried out here is only qualitative and to have 
a better understanding of the dynamics one needs to perform a thorough numerical calculation. 

One can also try to explore whether one can generate inflation [16] in these models without any 
explicit /9i nf i. Our preliminary numerical analysis does seem to suggest (see figure) that the external 
universe undergoes an accelerated phase of expansion at early times while the internal manifold is 
contracting. Such an inflationary regime will come to an end when branes annihilate to produce 
reheating. A related issue would be to study the stability of the branes in the first place. This 



t^Note that in our approach the equation of state for matter along the extra dimensions is the same as 
that of a pure cosmological constant term. 
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may reveal that a brane gas cannot be supported for some of the group manifolds and thus their 
smallness today cannot be explained in the BGC framework. 

Finally it would be interesting to try and compute the annihilation rates of branes wrapping dif- 
ferent subgroups (or coset spaces) of the internal group manifold and then apply the same reasoning 
as was done for the tori to see if one can generate a hierarchy of size (or squashing) within the full 
group manifold. It would be even more interesting to be able to start out with a non-compact simple 
group manifold and carry out this analysis to see whether the non-compact directions really become 
large while the compact ones are constrained by branes. This would then at least be able to explain 
why today we observe a product structure of the external and the internal universe. Of course this 
still will not explain the dimensionality of the observed universe. 
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